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ABSTRACT 



Context. Measures of the H I properties of a galaxy are among the most sensitive interaction diagnostic at our disposal. We report 
here on a study of HI profile asymmetries (e.g., lopsidedness) in a sample of some of the most isolated galaxies in the local Universe. 
This presents us with an excellent opportunity to quantify the range of intrinsic HI asymmetries in galaxies (i.e., those not induced 
by the environment) and provides us with a zero-point calibration for evaluating these measurements in less isolated samples. 

r ^ ' Aims. We aim to characterize the H I profile asymmetries in a sample of isolated galaxies and search for correlations between H I 

^^ I asymmetry and their environments, as well as their optical and far infrared (FIR) properties. 

^S^ . Methods. We use high signal-to-noise global HI profiles for galaxies in the AMIGA project (Analysis of the Interstellar Medium of 

QLh" Isolated GAlaxies, http : //amiga . iaa . es). We restrict our study io N - 166 galaxies (out of 312) with accurate measures of the 

HI shape properties. We quantify asymmetries using a flux ratio parameter. 

Results. The asymmetry parameter distribution of our isolated sample is well described by a Gaussian model. The width of the 
distribution is cr = 0.13, and could be even smaller (cr = 0.11) if instrumental errors are reduced. Only 2% of our carefully vetted 
^ ■ isolated galaxies sample show an asymmetry in excess of 3cr. By using this sample we minimize environmental effects as confirmed 

by the lack of correlation between H I asymmetry and tidal force (one-on-one interactions) and neighbor galaxy number density. On 
the other hand, field galaxy samples show wider distributions and deviate from a Gaussian curve. As a result we find higher asymmetry 
rates (~ 10-20%) in such samples. We find evidence that the spiral arm strength is inversely correlated with the HI asymmetry. We 
also find an excess of FIR luminous galaxies with larger HI asymmetries that may be spirals associated with hidden accretion events. 
Conclusions. Our sample presents the smallest fraction of asymmetric H I profiles compared with any other yet studied. The width of 

^^ I the associated asymmetry parameter distribution can help to distinguish the frequency and processes of self-induced HI asymmetries, 

\^ , and serve as a baseline for studying asymmetry rates in other environments. 
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1 ■ Introduction most extended cold component of the interstellar medium (ISM) 

and is a sensitive diagnostic of perturbations. The HI com- 

The geometry and kinematics of gaseous disks in galaxies ponent in spiral galaxies has long been kno wn to show both 

are mainly governed by the gravitational potential of stellar geometrical and kinematic asymmetries (e.g., "Beale & Davies| 

_ and non-baryonic components. Any perturbation of the equi- 1969. Ruchtm eier 1972 ; Allen et a l. 1973; B aldwin et a l. 1980). 

^ librium caused by either internal or external processes, may Asymmetries in stellar disks are also common and are traced by 

- - produce asymmetries m these disks. Atomic gas (HI) is the optical and near-infrared light (the latter less affected by dust 

extinction). Observations show that 30% of g alaxies are signifi- 

Send offprint requests to: D. Espada cantly lopsided at near-infrared wavelengths (Block et al." 1994 

* Based on observations with the 100-m telescope of the MPIfR iRix & Zaritskv 1995; Zaritsky & Rix 1997; Bournaud et al] 



(Max-Planck-Institut fuer Radioastronomie) at Effelsberg, GBT I2Q05I) . Asymmetries in the stellar component are not neces- 



under NRAO (the National Radio Astronomy Observatory is a sari ly correlated with lopsidedness in the gaseo us component 

facility of the National Science Foundation operated under co- (Ko rnreich et al.ll200Ql: IWilcots & Prescottll2004l) . This lack of 

operative agreement by Associated Universities, Inc.), Arecibo correlation is not surprising because the neutral hydrogen com- 

Observatory (National Astronomy and Ionosphere Center, which is p^^ent in a galax y is typically twice as e xtended as the stellar 

operated by Cornell University under a cooperative agreement with ^^ ^nt (e.g., iBroeils & van Wo"^M^ fT994h . and might be 

the National Science Foundation) and the Nancay Observatory, j ^ij-i-Sr . -.i ...i .n 

Full Table [D is available in electronic form at the CDS via anony- frturbed in diflFerent ways with respect to the stellar component 

mous ftp to cdsarc.u-strasbg.fr (13©. 79. 128. 5) or via (e.g., stripping). The study of HI properties is thus a better probe 

http : //cdsweb . u-strasbg . fr/cgi-bin/qcat? J/A+A/vvv/ppp of past and/or recent perturbations than the stellar counterpart, 

and from http : //amiga . iaa . es. especially for weak interactions. 
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Study of global H I velocity profiles of galaxies has proven 
to be very useful for a quantification of the frequency and ampli- 
tude of disk asymmetries (e.g. Richter & Sancisi 1994). While 
only aperture synthesis can provide full 2D information about 
the HI distribution and kinematics, the ID profiles provide a 
valuable measurement at a small fraction of the cost in observ- 
ing time. 

Past work suggests that the HI asymmetry properties 
of galaxies do not depend strongly on local environmen- 
tal conditions. Studies of field and/or isolated galaxies sug- 
gest t hat at least 50 % show significant HI profi le asymme- 
tries (iRichter & Sancisil 119941: iHavnes et aL 1998), an d eve n 
higher -75% in late-type spiral galaxies (IMatthews et al.lll998h . 
Although homogeneous studies of the asymmetry rate in richer 
environments are rare, it is usually believed that they show a 
comparable rate of asymmetric H I profiles. 

The implications of a high asymmetry rate essentially inde- 
pendent of e nvironment is ciirrently the subject of deb ate (for 
a review see iJog & CombesI l2009l: ISancisi et aP l2Q08h . It has 
been suggested that the mechanism responsible for producing 
asymmetric disks must be long-lived because high asymmetry 
rates are observed in samples of field and/or isolated galaxies. 
Because the signatures of tidal encounters are relatively short- 
lived, lasting only on the order of a dynamical time- scale for 
a wide range in mass ratios, orientations, inclinations, relative 
velocities and impact parameters (e.g. Bournaud et a_L 2005) , 
it cannot be the only agent responsible for the high asymme- 
try rate in diff'erent environments. A number of longer-lived 
mec hanisms have been proposed, a) intermi ttent minor merg- 
ers (^ Walkeretal1ll99 6': Zaritskv & Rix '199 7|), b) high- velocity 
cloud/gas accretio n (JBournaud et al. 2005: iSancisi et al.l 120081; 
Miller et aDl2008l). c) halo-disk m isalignment (JLevine & Sparkd 
19981: iNoordermeer et al.l l200ll) and/or d) internal perturba- 
tions including s ustained long-lived lo psidedness owing to non- 
circi ilar motions ([Baldwin et al.lll980l) or global m=l instabili- 
ties (ISaha et al.ll2007l) . 

In order to address the relevance of the diff'erent proposed 
mechanisms (internal versus external, short versus long-lived) 
one must first study a sample of well isolated galaxies in the 
nearby Universe (< 150 Mpc). This approach should minimize 
any contribution from tidal interactions and facilitate the inter- 
pretation of results with respect to other samples of galaxies. 
Reference samples used to study the rate of HI asymmetries 
involve galaxies which, although assumed to be field/isolated, 
usua lly include a significant population of interacting galaxies 
(e.g. iRichter & Sancisil [l 9941) . Note that field galaxies are de- 
fined as galaxies not belonging to the cluster environment and a 
significant number of th em are likely to be members of interact- 
ing pa irs or multiplets (ISulentic et al.ll2006b . IRichter & Sancisil 
(1 19941) find that about half of the nearby field galaxies show 
asymmetric profiles, estimated from a compilation of six H I sur- 
veys (1371 profiles were classified). It is important to remem- 
ber that, 1) the overall environmental properties of the samples 
were not assessed meaning that a significant fraction of galaxies 
might be environmentally influen ced, and 2) asymmetries were 
assessed using qualitative criteria (IRichter & Sancisi|[l994l) . 

Statistical studies of HI asymmetries in large samples of 
galaxies selected according to a well defined isolation criterion 
and using an objective quantification are rare. The only existing 
sy stematic study of H I asymmetries in an isolated sample is that 
of iHaynes et al.l (Il998l) . who studied the asymmetry rate for N = 
104 (A/^ = 78) galaxies that obey a0.5°(l°) projected separation 
criterion with respect to any known companion in the Arecibo 
General Catalog (AGC, private database of R. Giovanelli and M. 



R Haynes), where 0.5° corresponds to 175 kpc at a typical ve- 
locity of the core sample at Vr= 1500 km s"^ . About half of the 
galaxies appear to show significant H I asymmetries. However, 
only companions with a velocity diff'erence AV< 400 kms"^ 
relative to the primary are considered. The AGC is complete up 
to m ~ 15.4 mag and/or a diameter of V (equivalent to a typical 
linear size of 6 kpc). The velocity criterion could make this sam- 
ple biased against unbound plunging encounters, and owing to 
the size limit of the catalog, small galaxies would not have been 
taken into account - hierarchical systems of small galaxies may 
in principle produce significant asymmetries. 

Other studi e d sam ples are composed of field galaxies. 
IMatthews et al.l (Il998h studied a sample of A^ = 30 moder- 
ate to low surface brightness late-type spirals. Their galaxies 
lie between 2° and 6° from the center of the Fornax cluster 
or are field galaxies. They suggest that 77% of the HI pro- 
files in their sample show a relevant asymmetry. iBournaud et al.l 
(l2005h studied A^ = 76 galax ies based on the OSUBS Galaxy 
Survey (lEskridge et al.ll200"2l) . a sample not selected according 
to any environmental criterion, which hence may have a relevant 
amount of interacting galaxies. Howeve r, it has been recentl y 
referred to as a field galaxy sample (^e.g. Jjog & Combell2009l) . 
The asymmetry is larger than 10% for nearly 66% of the galaxies 
in this sample. 

In this pap er we use a large and complete sample of iso- 
lated galaxies (IVerdes-Montenegro et al.l 120051: iLisenfeld et al.l 
l2007l) to evaluate the intrinsic distribution of HI asymmetries. 
Comparisons show that our sample is more isolated than those 
used in previous studies and indeed is representative of the most 
isolated galaxies in the local Universe. The AMIGA project 
(Analysis of the interstellar M edium of Isolated GAlaxie£i 
IVerdes-Montenegro et al.ll2005l) involves vetting and analyzing 
the properties of galaxies i n the Catalogue of Iso lated Galaxies 
{N = 1050 galaxies, CIG. lKarachentsevalll973l) , and provides 
a good starting point for this aim. This project includes a re- 
fin ement of th e sample through 1) revision of optical positions 
(Le on & Verd es-Montenegro 20 ()3|), 2) analvsis of optica l prop- 
erties and completeness (Verde s-Montenegro et al.ll2005l) . 3) re- 



vised optical morphologies (iSulentic et al. 



uation of the isolation degree (IVerley et al.l 



2006h and 4) reeval- 



2007cllbh . We use afl 



these refinements in the present paper. The isolation criterion 
used for this compilation minimizes t he probability of a major 
interaction within the last ~ 3 Gyr (IVerdes-Montenegro et al.l 
l2005l) while quantifying possible minor interactions. A mul- 
tiwavelength characterization of diff'erent interstellar medium 
(ISM) components/phases and of the stellar component has been 
carried out including a) optical ( Verdes-Montenegro et al. 20051). 
b) FI R (ILisenfeld etal.ll2007 | ). c) radio-continuu m (iLeon et al.l 
l2008h . and d) U a emission (IVerleyetan 12007 ah . as weU as e) 
nuclear activity (ISabater et al.ll2008ll201^ . 



This paper presents an analysis of global (ID) asymmetry 
measures for N = 312 HI AMIGA galaxies with high signal-to- 
noise (S/N) spectra (Sect. [2]). We cleaned the sample of sources 
with uncertain asymmetry measures, yielding a total ofN= 166 
galaxies with high reliability data. For this subsample we studied 
the H I profile asymmetry rate (Sect. [3l), the role of the environ- 
ment for the rate of lopsided profiles (Sect. |4]), as well as the 
correlations between H I asymmetry and stellar properties (op- 
tical luminosity, morphological type and signs of perturbation), 
and star formation rate (as traced by FIR luminostity. Sect. [5]). 
Finally we identify the underlying distribution of intrinsic pro- 



http : //amiga . iaa . es 
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file asymmetries and discuss the possible origin of asymmetries 
in gaseous disks (Sect.©. 



2. The sample 

2.1. Sample selection 

This paper presents descriptions of the shape of HI profiles 
with special emphasis on the degree of asymmetry. We identi- 
fied the galaxies with the highest S/N H I profiles in the AMIGA 
project, obtained from observations at the Arecibo, Eff'elsberg, 
Nan^ay and GBT radio-telescopes, as well as from a compi- 
lation from archives/literature. The HI spectra are available at 
http://amiga.iaa.es. We initially consider for this study 
those H I spectra with a signal-to-noise S/N > 10 (S/N obtained 
as the peak flux to rms ratio), in total N = 3^3 galaxies. A S/N 
> 10 is appropriate to ga in a good estimate of the asymmetry 
(iTiff-t & Huchtmeiedll99d) . 

In addition we further restricted the sample to galaxies 
with profiles with high- velocity resolution to total width ratio, 
W20/AV > 10, where W20 is the width at a 20% level and AV 
is the resolution of the H I profile. This criterion is necessary 
to ensure that the profiles are well sampled. As a side-eff'ect, 
it partially excludes face-on and low mass/luminosity galaxies. 
We excluded sources with evidence of problems related to bad 
baseline subtraction and/or interference contamination. Finally 
we chose galaxies with recession velocities V > 1500 kms"^ to 
permit proper evaluation of the isolation. A prohibitively large 
region on the sky wou ld be required to asse ss isolation for galax- 
ies with V < 1500 (IVerlevetal.ll2007cllbh . In any case, they 
are all members of the Local Virgo Supercluster. Application of 
these restrictions yielded a sample of A/^ = 312 isolated galaxies 
with high quality H I profiles. We will refer to this sample as the 
HI sample. We characterize profile asymmetry in two ways: 1) 
quality-based on visual inspection of the profiles and 2) quantity- 
based on an areal asymmetry index Afiux ratio (Sect.O. For a sta- 
tistical analysis we performed a further refinement of the sample 
by considering only those galaxies with lowest uncertainty in 
their asymmetry parameter. This resulted in a final sample of N 
= 166 isolated galaxies, which we call the H I refined suhsample 
(Sect.[Il. 



2.2. Basic properties of the samples 

Figure [T] summarizes the basic optical/FIR properties of both 
the HI sample and the HI refined subsample. The diff'erent 
panels include recession velocity, morphological type as well 
as optical and FIR luminosity distributions. We compare these 
distributions with those corresponding to the optically com- 
plete sample, which was estimated to be 85-90 % complete 
to mB= 15.0 mag, and is comp osed of N = 719 GIG galaxies 
(IVerdes-Montenegro et al.ll2005l) . The latter sample is appropri- 
ate to represent the entire population of isolated galaxies in the 
local Universe. 

We summarize the dispersions of the basic properties for 
galaxies included in the diff'erent samples as well as the devi- 
ations between the H I samples and the optically complete sam- 
ple: 



skewed toward lower velocities relative to the optically com- 
plete sample. The sample becomes seriously incomplete be- 
yond 9000 km s" 1 . 

- b) Morphology T(RC3) (ISulentic et al.l l2006h: Types are 
given in the RG3 numerical scale (Ide Vaucouleurs et aP 
I1991I) . The bulk of the GIG sample involves late-type galax- 
ies in the range 3 <T <1 (Sb to Sd) with 2/3 of the sample 
in a very narrow range T=4±l (Sb-Sc). Only 14% of the 
sample involve early-type systems, which suggests that our 
sample represents the extreme (low) end of the morphology- 
density relation. The H I refined subsample contains a higher 
percentage of late-type galaxies (especially Sb-Sc) than the 
optically complete sample. This is not surprising because 
earlier type galaxies have a systematically lower HI con- 
tent than later types and are especially excluded by selecting 
those H I profiles with S/N > 10. 

- c) Optical luminosity Lb (IVerdes-Montenegro et al.ll2QQ5h : 
With few exceptions the sample spans a 1 dex lumi- 
nosity range (9.5 < log(L5[Lo]) < 10.5). Galaxies with 
log(Lfi[Lo])< 10 are overrepresented in the HI sample 
(and HI refined subsample) likely because higher lumi- 
nosity galaxies prefer the high-velocity tail of the sample 
(Malmquist eff'ect) and often fall below m= 1 5.0. 

- d) FIR luminosity Lpm (iLisenfeld et al.l20()7h : The three dis- 
tributions are relatively similar to each other. The shape of 
the FIR luminosity distribution is flatter than the optical and 
shows a peak near log \^fir= 9.6. The full range covers 2 dex 
(8.5 < log(L/7//?[L0]) < 10.5). The difl'erence between the op- 
tical and FIR luminosity distributions likely reflects the ex- 
treme FIR "quietness" of our very isolated galaxy sample. 
Only detections are shown in Figure [T] 



3. H I profile shape and quantification of H I 
lopsidedness 

In this section we present a general view of the H I profile shape 
and two ways of quantifying profile lopsidedness. First, we ex- 
amine the profiles via visual inspection (see Sect. 13.11) using cri- 
teria similar t o those employed in the l argest H I profile shape 
study to date, iRichter & Sancisil (Il994l) . We then quantify the 
asymmetry level in a more objective manner using a numerical 
parameter (Sect. 13.21) and compare visual and numerical descrip- 
tions in Sect. 13.31 



3.1. HI profile stiape and visual estimation of lopsidedness 

Visual inspection of our HI sample shows that 88% show 
double horns with the rest showing single peaks that in most 
cases involve face-on spiral galaxies. We visually classified the 
profiles in the HI sample (N = 312 galaxies) as symmetric, 
sl ightly asymmetric, and st rongly asymmetric, in a similar way 
as 



Richter & Sancisil (11994 . who studied a sample of N = 1371 



- a) Radial velocity V (IVerdes-Montenegro et al.l l2005h : 
Velocities range over Idex in the three samples, covering the 
range 1500 kms"^ < V < 14000 kms"^. The HI refined 
subsample, and to a lower extent the H I sample, is slightly 



spectra (equivalent, respectively, to their notation as No, Weak 
and S trong). We find that A^ = 141 galaxies show symmetric H I 
profiles (45 %), A/^ = 126 sligthly asymmetric profiles (40%), and 
N = A5 strongly asymmetric profiles (15%). In order to illustrate 
this visual classification we show some examples of symmet- 
ric, slightly asymmetric, and strongly asymmetric H I profiles in 
Figures [2I O andlH respectively. Based on the visual classifica- 
tion, our sa mple appears to show sim ilar average rates as those 
obtained bv lRichter& Sancisil (11994 : 47 + 5 %, 34 + 6 %, and 
19 + 6 %. 
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The most common asymmetry found in our sample involves 
unequal peaks in the double horn profiles. The most extreme 
HI asymmetries occur for 43 galaxies (out of 312) where the 
peak flux difl'erence is approximately larger than 25%. Only one 
of them has a peak flux difl'erence larger than 50% (CIG 317). 
Some galaxies with double peaked profiles show peculiarities: 
CIG 144 shows a central peak stronger than the horns of the dou- 
ble peaked profile, CIG 858 shows a profile with a peculiar cen- 
tral trough, CIGs 238, 382, 928, and 1029 show apparent wings 
(significant excess flux -3cr- within a 50 - 100 km s"^ wide) be- 
yond the outer walls of the double-horn profile. CIG 170 shows 
an uncommon flat H I profile. CIG 870 may also have wings that 
are 50 - 100 kms"^ wide, although it seems to be a face-on 
galaxy. The observed wings may indicate a projected gas-rich 
companion or extra-planar motions owing to a nurture event. 



3.2. Integrated density flux ratio parameter (Afiu, 



o) 



A variety of parameters have been used in the literature to quan- 
tify the asymmetry level in HI profiles. We employ an areal 
asymmetry index to quantify profile lopsi dedness, namely the 
integrated flux densit y ratio Afiux ratio (e.g. iHaynes et al.lll998[ 
iKornreich et al.ll200ll) . defined as Afiux ratio = Mih^ '^^ Mih > 1, 
and 1 lAiih otherwise, where Aijh is the ratio of the areas under 
the profile at velocities lower (S/) and higher (S/^) than the central 
velocity (Vm)'- 



A - Si _ yi 



Vm 



Sh 



fS,,dv 



where V/ and Vh represent the low and high velocities measured 
at 20% intensity level with respect to the peak. Vm is calculated 
as the mean velocity at the same level, Vm = (Vh+Vi)/2. Note that 
^fiux ratio is invariant to the sense of rotation of the galaxy. We 
use this parameter since it is the most common asymmetry index 
that can be found in the literature and allows us to compare our 
results with other samples of galaxies (see Sect. 14.21) . Equivalent 
definitions are found in the bibliography and can be easily con- 
verted to Afiux ratio, as e.g.! A = ^^ (iMatthews et alJil998) and 
^1 = 10 X (1 - l/Afiux ratio) (iBournaud et al.ll2005h . We indi- 
cate the Afiux ratio valucs in Eigure [2l [3l and |4] for the examples 
of H I profiles visually classified as symmetric, CIG 226, slightly 
asymmetric, CIG 421, and strongly asymmetric, CIG361, which 
are characterized by Afiux ratio = 1-05, 1.15, and 1.51, respec- 
tively. 



3.2.1 . Uncertainties of the Afiux ratio 

We estimate the uncertainty of the asymmetry index, AAfiux ratio, 
by taking into account a) the rms noise per channel, b) the un- 
certainty in the calculation of the mean velocity, and c) the ob- 
servational /?6>m tog offsets: 

- a) Uncertainty owing to the rms of the H I profile (AA(rms)): 
owing to the rms of the spectrum, cr, the uncertainty in 5" / is 
A^" / = ^[NicrR, where Ni is the number of channels corre- 
sponding to Si and R is the spectral resolution of the profile. 
The uncertainty mSu, ASu can be obtained in the same way. 
Then AAfiux ratio Can be calculated as 

^flux ratio = KJ^^Sif + (|A^,)2)]1/2. 

- b) Uncertainty owing to the measurement of the mean ve- 
locity (AA(mean vel)): since the Afiux ratio is calculated as an 
areal ratio obtained from the mean velocity, an error in the 



determination of the latter can induce a wrong measure of 
^fiux ratio- An crror AVm in the estimate of the mean velocity 
Vm produced by limited velocity resolution and/or S /N ratio 
would artificially increase the asymmetry index of a symmet- 
ric profile (i.e. Afiux ratio = 1). Namely Afiux ratio = f^ > 1. 

A good estimate of the uncertainty can be obtained from the 
increase in Afiux ratio for a symmetric profile. If the profile 
were symmetric, then S/ = S/^ = S/2, where A is the total area 
under the profile, e can be estimated as 6 ~ /z AVm, where 
h is an intensity height scale. The uncertainty of AVm can 

be estimated as AVm = 4-^^ (JEouque et al.lll99Cih , where 
P = (W2o-W5o)/2, parameter that represents the steepness of 
the edges of the H I profile, and W20 and W50 are the widths at 
20% and 50% with respect to the peak, respectively. Owing 
to the uncertainty in the determination of V^, we can express 

AS I (and ASh) as: AA/ = AAh = 4-^^/z, where we estimated 
that h = hmaxf^, being hmax the H I profile strongest peak. 
- c) Uncertainty owing to pointing off'sets (AA(pointing off- 
set)): a pointing off'set of the antenna with respect to the 
kinematic center of the galaxy can induce an artificial lop- 
sidedness in the H I profile when the tele scope beam is com- 
parable to the size of the galaxy (e.g. iTifl't & Huchtmeie j 
ll990l:ISpringob et ani2005l) . In some cases the observing co- 
ordinates were not coincident with the center of the galax- 
ies owing to errors in the position s found in the CIG (e.g. 
iLeon & Verdes-Montenegroll2003l) . The expected flux loss 
(f) owing to beam attenuation and known antenna pointing 
offsets can be calculated from t he optical dianieter o f the 
galaxy, beam size, and pointing (ISpringob et al.ll2005h . We 
decomposed the expected flux loss (f) into two components, 
/ = fb.a. fp.o. , where /^.«. is the flux loss factor arising from 
beam attenuation and fp_o. is the factor where the contribu- 
tion of the pointing ofl'set is. The latter is intimately related 
to the flux loss that contributes to the asymmetry of the H I 
profile. Note that beam attenuation with no pointing oflfset 
causes no asymmetry in the H I profile. The diff'erence in the 
^fiux ratio parameter if the flux loss that contributes to the 
asymmetry is on the receding or approaching side provides a 
measure of AA(pointing offset). We find among the HI pro- 
files only eight galaxies (out of the 312) that have diff'erences 
larger than 0.01. The average is a factor 10 smaller than the 
contribution from the other two sources of uncertainty. The 
average is 0.001 and the standard deviation is 0.004. Thus, 
this source of error is negligible in most cases in our data. 
This is a result of the small known pointing off'sets in the ob- 
servations, where the average is equal to 3'' and the standard 
deviation is 3'\ 

We added these sources of uncertainty in quadrature to esti- 
mate the net uncertainty in Afiux ratio, ^^fiux ratio- 

3.2.2. Other possible sources of uncertainty 

In addition to these sources of uncertainty, there are other eff'ects 
that can induce an artificial asymmetry on the H I profiles. These 
include the eff'ect of random pointing off'sets and baseline fitting. 
Elux l oss due to this efi^ec t must be somewhere between < 
1% (GBT. iHavnes et al.lll998l) and 5% (Arecibo circular feed, 
iHaynes & Giovanellil 19841) . If the profile were initially symmet- 
ric, then the induced asymmetry parameter by this eff'ect would 
be in the range Afiux ratio < 1.02 - 1.1 1, if all flux loss is located 
either in the receding or approaching sides. Because we have 
data from difl'erent telescopes, our situation is probably interme- 
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diate between both cases. We assume that the resulting Afiux ratio 
distribution of a sample with symmetric H I profiles observed 
under similar conditions as our sample is likely well represented 
by a half-Gaussian with a cr = 0.04. 

The baseline fitting proc ess can also produce a rtificial asym- 
metrie s in the HI profiles (iHaynes et al.l Il998b . iHaynes et al.l 
(1 19981) indicate that diff'erent order fits show flux diff'erences of 
about 3%. As a result, the asymmetry parameter for symmet- 
ric HI profiles can be altered up to Afi^x ratio = 1-06. A half- 
Gaussian curve with cr = 0.02 would mimic this eff'ect well. 

Given the random nature of these two eff'ects, we cannot es- 
timate their values individually, but their overall eff'ect is taken 
into account in Sect.[6lto discuss the actual shape of the Afiux ratio 
distribution in isolated galaxies because they can broaden the re- 
sulting distribution. 

3.2.3. Presentation of the asymmetry data 
We list in Table [T] the following information: 

- 7j GIG number; 

- 2) Visual classification: = symmetric, 1 = slightly asym- 
metric, and 2 = strongly asymmetric (Sect. l3.ll) , 

- 3) Afiux ratio, the asymmetry parameter (Sect. l3.2l) , 

- 4) AA(rms), the uncertainty in Afiux ratio owing to the rms of 
the HI profile (Sect. [ITB, 

- 5) AA(mean vel), th e unce rtainty in Afiux ratio owing to the 
mean velocity (Sect. l3.2."T]) , and 

- 6) AA flux ratio, the global uncertainty in Afiux ratio, in- 
cludi ng the small contribution from AA(pointmg offset) 

The Afiux ratio distribution is shown in Figure \5\ The best 
half-Gaussian fiQ to the asymmetry parameter distribution is 
characterized by a cr = 0.15. However, this half-Gaussian fit is 
not able to reproduce the Afiux ratio distribution both at the high 
and low ends. First there is an excess of high values of Afiux ratio 
with respect to the Gaussian curve, and second, the peak of the 
distribution is too flat for Afiux ratio < 1-15. 

We show in Figure [6l a and b the AA(rms) and AA(mean vel) 
distributions, respectively. The combined efl'ect of all the previ- 
ous uncertainties, AAfiux ratio (including the small contribution 
of AA(pointing offset)), is shown in Figure [6] c. We show the 
best Gaussian fits to the distributions. 



3.3. Comparison between the visual classification and 

^flux ratio 

We compare the asymmetry visual classification of the H I pro- 
files (Sect. 13.11) with the Afiux ratio in Figure [71 Three clearly 
distinct Afiux ratio distributions are seen for those galaxies visu- 
ally classified as symmetric, slightly asymmetric, and strongly 
asymmetric (Sect. 13.11) . The Afiux ratio distribution of HI pro- 
files visually classified as symmetric has a mean value equal to 
1.08, with a standard deviation of 0.065. The distribution for the 
slightly asymmetric H I profiles is characterized by a larger mean 
of 1 . 13 and a similar standard deviation, 0.09. The distribution of 
strongly asymmetric profiles is characterized by a mean of 1.37 
and a considerably larger scatter, 0.17, with values as high as 
Afiux ratio = 1.8. The Afiux ratio distribution for the sHghtly asym- 
metric subsample partially overlaps with those of the symmetric 
and asymmetric distributions. 



^ We fitted the parameters A, fi and cr in a half-Gaussian curve defined 

as|Aexp(^^)| 



The large overlap that exhibits the Afiux ratio distribution for 
each visually classified subsample is not surprising, because this 
visual classification is of course subjective, and because the 
Afiux ratio parameter misses a few cases where the shape of a 
real asymmetric profile does not correspond to diff'erent areas in 
the approaching and receding sides. Future work would require 
to inspect other asymmetry parameters that are sensitive to flag 
these profiles as asymmetric. 



3.4. HI refined subsample and characterization of the 
Afiux ratio distribution in a sample of isolated galaxies 

The shape of the Afiux ratio distribution might be aff'ected by arti- 
ficially induced values from the eff'ects explained in Sect. 13.2.11 
By reducing the net uncertainty in the asymmetry measurement, 
we reduce errors that might bias our results. We show in Figure[8] 
how the Afiux ratio distribution changes for diff'erent AAfiux ratio 
limits. The smaller the limit (i.e., only including accurate values 
of Afiux ratio), the better a half-Gaussian reproduces the distribu- 
tion. 

From now on we choose those H I profiles with AAfiux ratio 
< 0.05, namely the HI refined subsample, to remove from our 
statistical analysis those profiles with an uncertain determina- 
tion of the asymmetry index. With this criterion we still have 
a large sample of N = 166 galaxies. The basic property distri- 
butions (velocity, morphological type. Lb and Lfm) of the HI 
refined subsample are shown in Figure[T]as (blue) solid lines, in 
comparison to those of the H I sample. 

In order to characterize the intrinsic scatter of the asymmetry 
parameter distribution in a sample of isolated galaxies with mi- 
nor contamination of artificially asymmetric H I profiles we fit- 
ted a half-Gaussian function to the H I refined subsample (FigO. 
The fit yields a width of cr = 0.13 (Fig [9]). This time the fit suc- 
cessfully reproduces the asymmetry parameter distribution, in- 
cluding the low and high ends. Only 2% of the isolated galaxies 
are in excess of 3(T. 

The width of the half-Gaussian distribution sets an upper 
limit to the intrinsic dispersion of the H I asymmetry in isolated 
galaxies. Errors in the calculation of the asymmetry index might 
be typically ~ 0.03 (mean of the Gaussian fit) (Sect. 13.231) . As 
discussed in Sect. l3.2.2l there might also be random errors in the 
pointing (~ 0.04) and baseline subtraction (cr ~ 0.02) that may 
increase errors in Afiux ratio- Hence it is reasonable to expect a 
lower value of the width, cr ~ 0.1 1, once these sources of errors 
are corrected. 

Note that the quantification of the asymmetry distribution for 
the galaxies in the H I refined subsample is not aff'ected by incli- 
nation eff'ects (e.g. I Jog & Combesll2009l) . Figure [TOb shows that 
the inclination of the galaxies are distributed homogeneously 
above / = 30°. Only two galaxies have an inclination i < 15° 
(GIG 85 and 178). The lack of galaxies below / = 30 is caused by 
the width-to-channel ratio criterion explained in Sect. 12.11 This 
homogeneity in the inclination ensures that most of our galaxies 
do not show symmetric profiles because the galaxies are viewed 
face-on, where an asymmetry in the velocity field would remain 
unnoticed. To further inspect whether inclination can be intro- 
ducing any bias in our results, we plotted the asymmetry index 
versus the inclination (Figure [TOb) and found that the two quan- 
tities are not correlated. 
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4. H I profile lopsidedness and environment 

4.1. HI asymmetries and isolation parameters in CIG 
galaxies 

A reevaluation and quanti fication of isolation deg ree for CIG 
galaxie s was reported in IVerley et aP (l2007cllbb . IVerley et aP 
(l2QQ7d) derived two isolation parameters for each CIG galaxy: 
1) a local surface density parameter t]k within the distance to the 
k-th neighbor (a good tracer of average galaxy surface density) 
and 2) a tidal strength parameter Q (a parameter more sensitive 
to one-on-one interactions). 

H I is known to be a sensitive diagnostic of interaction mo- 
tivating us to compare these two parameters with our Afi^x ratio 
asymmetry parameter. Figure \TT\ shows the lack of correlation 
between Afiux ratio and both t]k and Q. The Pearson's correla- 
tion coefficient is p = -0.005 and 0.114, respectively, which in- 
dicates that the two quantities are essentially not correlated. A 
small trend in the Q parameter might be present in the sense that 
larger H I asymmetries seem to be found in less isolated systems. 
The calculated intercept and slope are -3.3 ±0.7 and 0.8 ± 0.7, 
respectively. 

The lack of correlation suggest that we are minimizing nur- 
ture eff'ects that might aff'ect the H I shape. The low values and 
small range in terms of galaxy density and tidal strength covered 
by CIG galaxies are not enough to see a correlation. 

iBournaud et al.l (l2005h also suggest that there is no correla- 
tion between lopsidedness and tidal strength. However, they use 
the lopsidedness Ai parameter on NIR surface density, and NIR 
emission is not as extended as the H I. 



4.2. H I asymmetry distribution in field samples 

We compare the asymmetry distribution of our H I refined sub- 
sample with that of diff'erent studies from the bibliography where 
a similar asymmetry index (Sect. 13.21) has been c alculated and 
involv e field/isolated galaxies (see also Sect. [B : Havnes et al.l 
(ll998h . lMatthews et al.l (Il998h and lBournaud et al.l (l2005h . 

Figure[T2l shows the A fjux ratio normalized dist ribution for our 
refined subsample with D iBournaudet al.l (|2005|) and 2) a com- 
bined s ample (N = 1 86) incl udin g H I data in iMatthews et al.l 
(ll998l) . lBournaudetal.l (l2005h . and lHavnes et al.l (Il998l) exclud- 
ing CIG galaxies (80 galaxies). The H I refined subsample shows 
the distribution best described by a half-Gaussian It also shows 
the lowest absolute value of cr. The lBournaud et al.] (l2005l) distri- 
bution shows the widest distribution (cr = 0.23) and noticeably 
deviates from a half-Gaussian curve. An intermediate case, cr 
= 0.17, is found for the combined sample without CIG galax- 
ies. Table [3] gives cr values for each distribution as well as an 
asymmetry rate with "asymmetric" profiles defined as Apux ratio 
values exceeding the 2cr level of our HI refined subsample 

{A flux ratio — 1.26). 

Figure[T3]compares the Afi^x ratio cumulative p robability dis- 
tributi o n for our HI refined sa mple and those of iHaynes et al.l 
(119981) . IMatthews et aP (Il998h . and iBournaud et al.l (l2005h . In 
each plot the diff'erence of the two curves indicates the asym- 
metry rate diff'erence for a given Afiux ratio Hmit. Our sample lies 
below the field samples in almost every bin with diff'erences typ- 
ically bet ween 10 - 20%. A re sult more similar to our sample is 
found for lHaynes et al.l ([1998') likely in part because of the sig- 
nificant fraction of CIG galaxies (23%) included in their sample. 
Removing the CIG overlap increases both their cr and asymme- 
try rate. 



We performed a x^ test to check whether the null hypothe- 
sis that any of the three Afiux ratio distributions is similar t o our 
H I refi ned sample, could be rejected. Except for Haynes et al.l 
(Il998h (x^ = 9 and the associated p -value = 0.33) this hyp oth- 
esis can be reje cted. In the cases of IBournaud et al.l (l2005l) and 
IMatthews etakl ([1998), x^ = 47 (p-value = 2 x 10"^) and x^ 
= 14 (p-value = 0.09) respectively. The sample diff'erences we 
find are significant and cannot be ascribed to the refinement of 
the H I sample (Sect. 1341). I n pri nciple we do not know how 
much IMatthews etaP (Il998h and IBournaud eFaP (l2005h 's ob- 
servations are aff'ected by systematic errors, but diff'erences in- 
volving the same criterion as use d in our study would yield an 
asymmetry rate diff'erence < 5%. iHaynes et all (Il998l) included 
only high S/N profiles, avoided pointing problems and quanti- 
fied baseline problems, suggesting it is reasonable to compare it 
directly with our H I refined sample. Overall, bec ause of their de - 
gree of isolation, our HI refined subsample and IHaynes et al.l 's 
show a lower frequency (~ 10 - 20%) of galaxies with asvmmet- 
ric profiles than in other sa mples such as IMatthews et al.l (Il998l) 
and lBournaud et al.l (l2005l) . 



5. Relation between H I profile lopsidedness and 
optical/FIR properties 

In this section we explore possible correlations between the 
asymmetry index Afiux ratio and optical properties of the H I re- 
fined subsample such as morphology, optical signs of perturba- 
tion, optical luminosity (Lb), and far-infrared luminosity (Lfir). 

5.1. Morphology and luminosity 

Figure [14] shows the distribution of Afiux ratio values for each 
Hubble morphological class (median, mean, and standard devia- 
tion values are indicated). The bins representing the majority of 
our sample (T(RC3) = 3 to 6, i.e. Sb to Scd) show a fairly large 
scatter (standard deviation cr^ 0.1). We see a slight decreasing 
trend in Afiux ratio toward later- type galaxies. 

Studies of the relati on between H I lopside dness and mor- 
phological type are rare. IMatthews et al.l (1 19981) studied a sam- 
ple of 30 moderate to low surface brightness (T=6-9) galax- 
ies and found a higher asymmetry rate than for more luminous 
(and higher surface brightness) late-type spirals. They found (for 
this type range) that later types were more likely to show larger 
asymmetries. In the Eridanus group the Ai parameter as cal- 
culated from H I maps is larger for earlier type galaxies, sug- 
gesting that tidal interactions generate a higher lopsidedness 
rate in galaxies undergoing secular evolution toward earlier type 
( Angiras et al. 2006). This result (i.e., larger asymmetries for 
earlier types) agrees with the general trend seen in our sample. 

We compared the Afiux ratio with luminosity. More lumi- 
nous galaxies are slightly more asymmetric (Figure \T5\ left 
panel). Figure [15] (right panel) presents the cumulative distribu- 
tion of Afiux ratio fo^ the high- a nd low-luminosity subsamples 
(IVerdes-Montenegro et al.|[2005h . The two distributions are dif- 
ferent at a level of a = 0.05 using a chi square test: ;^^=14 and 
p - value=O.Ol. 

5.2. Optical signs of interactions 

Here we inspect a possible connection between optic al signs of 
interaction and asymmetries in the HI profiles. Sulentic et al.l 
(2006) revised the optical morphology classification for the 
CIG sample using P0SS2/SDSS data. Although the CIG, the 
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starting sample of AMIGA, has been selected to minimize 
clos e neighbors to the target galaxy and thus interactions, still 
ISule ntic et aL/s revision revealed N = 193 objects with nearby 



measures, including upper limits. Figure [T7] (left panel) presents 
a slight trend in the sense that more luminous FIR galaxies 
have more asymmetric HI profiles. Right panel of Figure [TTl 



com panions or signs of d istortion likely caused by an interac- shows the cumulative probability distribution for Afiu, 



com panions or signs oi d istortion liKely caused Dy an interac- 
tion. Hulenticelal] (l2006h flagged these galaxies as interacting 
in the case of a morphologically distorted system and/or almost 
certain interacting system or flagged as possibly interacting if 
there was any evidence of interaction/asymmetry with/without 
certain detection of a close companion. 

There is no statistically significant diff'erence in terms of H I 
profile asymmetry rate between galaxies that are optically clas- 
sified as interacting and those without any sig n of interaction 
This r esult is consistent with the conc lusion by iKornreich et al.l 
(I2QQQI) and lWilcots & PrescottI (|2Q04 that optical asymmetries 
are not necessarily correlated with a lopsided H I component. 



5.3. Bar and spiral strengths 

We compared the H I asymmetry parameter with the relative spi- 
ral and bar strengths calculated as the maximal torque, or ratio 
of the maximum tangential force and the azimuthally averaged 
radial force. This h as been obtained for a subsample of N = 
96 GIG galaxies by iDurbala et al.l (l2009l) usi ng Fourier analy - 
sis over spiral and bar components separately (JButa et al.ll20Q3l) . 
Figure [16] shows the Afiux ratio parameter with respect to the rel- 
ative spiral strength (Qs)- The overlapping sample between our 
H I refined subsample and the one used by Durbala et al. (2009) 
is composed of 40 galaxies. Qs seems to anti-correlate with 
A flux ratio ' disks with weaker spiral arms show stronger asym- 
metries (Figure O. There are six galaxies (GIGs 11, 33, 689, 
712, 912, and 931) that are outliers to this relation in the high 
end of the Afiux ratio parameter. These are likely galaxies whose 
H I asymmetry parameter is aff'ected by instrumental eff'ects. On 
the other hand we do not find any correlation of Apux ratio and 
Qb. The total strength Qg is not correlated with Afiux ratio ei- 
ther, not surprisin gly because Qh and Qg presents a good correla- 
tion (iDurbala et al. 2009) . The relation between Qs and Afiux ratio 
may originate in the observed trends in Qs versus T(RC3), be- 
cause Qs would be expected to correlate with the latter. However, 
we did not find a clear trend between Qs and T(RC3). 

The correlation found between Afiux ratio and Qs suggests 
that other samples of galaxies characterized by lower spiral 
strengths will have higher HI asymmetry rates. T his is consis- 
tent with the OSU sample (lEskridge e t al.ll20()2b , where spiral 
strength is lower tha n for the GIG (median equal to 0-132 v ersus 
0.161, respectively) (iButa et al.ll2"005l: IDurbala et al.ll2009l) . and 
have a higher H I asymmetry rate with respect to the H I refined 
subsample (see Sect. 14.21 and iBournaud et al.ll2005h . 



5.4. FIR luminosities (Lpir) 

LfjR is a good tracer of the star-formation rate and is re- 
lated to the environment in the sense that IR luminous galax- 
ies (LfiR > 10^^ Lq) are u sually interacting or merger systems 
(ISanders & Mirabel! 1 19961) . Unlike other samples of galaxies, 
our isolated population shows low FIR measures e.g. log(LFiR) 

S )eaks from 9.0 - 10.5 with very few (<2%) galaxies above 10.5 
Lisenfeld et al.l 120071) . The low LfjR values of the GIG sam- 
ple support our claim that the r evised GIG (AMIGA) is a sam- 
ple with only isolated systems (JLisenfeld et al.ll2007h . Here we 
inspect whether the small fraction of IR-luminous systems in 
our sample corresponds to galaxies with larger HI asymmetries. 
We find N=165 galaxies in our refined HI sample with IRAS 



three bins: 8.0 < log(L/.//?[Lo]) < 9.5, 9.5 < logiLpmiLQ]) < 10.0 
and 10.0 < log(L/7//?[Lo]) < 11.0. We also show the distribu- 
tion for those galaxies (N= 60) with FIR upper limits. The 
null-hypothesis that the first two bins are similar cannot be re- 
jected (at a level of a = 0.05) using a;^^ square test: ;^^=6 and 
p - value=0.24. On the other hand we find that the Apux ratio 
distributions for the latter two bins are diff'erent (x^=\A, p - 
va/w^=0.008) from each other. We find a 10-20% excess of 
higher asymmetry values for the most FIR luminous galaxies 
(10.0 < log(L/7//?[Lo]) < 11.0). If real, this excess might reflect 
asymmetries and FIR luminosities simultaneously enhanced by 
accretion events. 



6. Discussion 

The GIG sample represents the -3% most isolated galaxies 
in the Gatalqg of Galaxies and Glusters of Galaxies (GGGG, 
IZwicky et al.ll 19611) . In these systems the efl'ects of environmen- 
tal perturbation such as tidal interactions and ram pressure strip- 
ping are minimized. The Afiux ratio distribution of our sample 
is well described by a half-Gaussian function with a Icr width 
of 0.13 (possibly smaller cr ~ 0.11 if artificially induced eff'ects 
are reduced). Only 9% of the galaxies show Afiux ratio > 1.26 
(2(t) and only 5% Afiux ratio > 1-39 (3cr). If it is reasonable to as- 
sume that the distribution of intrinsic asymmetries should show a 
Gaussian distribution, then that of Figure[9]is as close as we have 
ever come to isolating that intrinsic distribution. Gomparison 
with field samples clearly shows that eff'ects of nurture result 
in an excess population of high Afiux ratio values. This excess 
population measured as a deviation from the best Gauss fit is 
negligible in our sample. In isolation it is apparently very un- 
likely to find galaxies with HI disks showing Afiux ratio > 1-39. 
The small number of such extreme asymmetric profiles found 
in our sample show double-peaked profiles with unequal horns. 
This is in most cases not caused by contamination by gas-rich 
companions with systemic velocity suitable to create a false or 
amplified horn. We almost never observe such components in the 
middle or close to the edges (thus broadening one of the horns 
in) the H I profile. We are unlikely to find a narrower asymmetry 
distribution in any galaxy sample. 

Unfortunately, using HI profiles does not allow us to dis- 
tinguish between the roles of geometry and kinematics in pro- 
ducing an asymmetry. Aperture synthesis maps suggest that 
asymmetry is usually the signature of kinematic lopsidedness 
although galaxies with lopsided HI distributions are not un- 
known. The rotation curve on one side of the ga l axy is usu- 
ally steeper than the opposite o ne. IS waters et al.l (1 19991) esti- 
mat ed from previous HI maps (Broeils & van Woerdenlll994l: 
iRhee & van Albadal [19961: IVerheiien 1997h that the fraction of 
kinematically lopsided galaxies may be as large as 15-50%. It 
is likely that the few asymmetric H I disks in our sample present 
this kinematically lopsidedness. However, to make a quantita- 
tive relation, it would be necessary to calibrate (statistically) how 
asymmetry parameters in 2D maps relate to those ID parameters 
using single-dish data. 

We searched for correlations between internal properties of 
galaxies and the measured Afiux ratio parameter. Although we 
maximized our sensitivity to internal correlations by removing 
all galaxies likely to have been aff'ected by external perturbers. 
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we did not find any strong correlations between stellar proper- 
ties and HI asymmetry. We find a weak correlation between 
spiral arm strength and H I asymmetry parameter, in the sense 
that arms are stronger for galaxies with more symmetric H I pro- 
files. The simplest interpretation for such a trend would be that 
gas asymmetries are more efficiently suppressed by the stronger 
spiral arm gravitational torques in more massive galaxies. A 
connection to spiral arm strength rather than to bar features 
should arise from the larger scale of the former. It has long been 
known that more developed spiral arms seem to exist in more 
luminous galaxies (Ivan den Berghlll96Ql) . Interestingly, the mor- 
phology of spiral arms is found to depend prim arily on parent 
galax y properties rather than on the environment (Ivan den Berghl 
l2002h . Therefore we would expect an anti-correlation between 
H I asymmetry and luminosity. However, Fig.[T5]shows a hint of 
the opposite trend. This might be because unlike van den Bergh's 
(1960) study, our study focuses on a sample with a small range 
in luminosity, as shown in Figure [T]:. 

The lack of a strong correlation between L^ir and Afiux ratio 
indicates that the bulk of the star formation and the symmetry 
of the gaseous disk are not strongly linked, i.e., that induced 
SF caused by possible interactions in this sample is small com- 
pared to that from secular evolution. Still, there is an excess 
of about 10% of asymmetric profiles for the most IR luminous 
(10 < log(L/7//?[Lo]) < 11) galaxies. This might be linked to re- 
cent accretion events in a small number of CIG galaxies. 

In general the asymmetry distribution will likely deviate 
from a half-Gaussian curve for other samples containing galax- 
ies that are perturbed by the environment. These samples are the 
rule while very isolated galaxies are the exception. The intrinsic 
asymmetry distribution found in our sample of isolated galaxies 
will be skewed toward higher values as a result of these inter- 
actions. This is confirmed by the wider distributions found in 
samples of field galaxies (Sect. 14.2b where a higher degree of 
interaction is expected, given the lack of a strict isolation cri- 
terion. Tlie_de}n^ti£nfrom_a half-Gaussian curve for the sam- 
ple in iBournaud et al.l (l2005l) is apparent as shown in Sect. 14.21 
and the distribution is the widest (cr = 0.23) among those stud- 
ied. Although it has been known for a long time that interacting 
galax ies usually show larger HI asymmetries (ISulentic & ArpI 
Il983l) . a statistical analysis using a common H I profile asymme- 
try parameter in large and well-characterized samples in dense 
environmen ts is s ti ll ne eded. 

Mapelli et al.l (1200 8) estimate from the density of (colli- 
sional) ring galaxies in the local Universe that major asymme- 
tries for ~ 10% of the galaxies may be produced as a result 
of a recent fly-by, resulting in a lopsidedness visible over a 
time scale of 1 Gyr. Within the uncertainties of this estimate, 
this may well match the 10-20% difl'erence between our sam- 
ple and field samples. As reviewed in Sect. [T] many mecha- 
nisms have been proposed other than tidal interactions or ram 
pressure from the intergalactic medium that may contribute to 
the asymmetry par ameter distribution, su ch as minor interac- 
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internal perturbations including sustai ned long-lived lopsid ed- 
ness owing to non-circ ular motions (JB aid win et all Il980l) or 
global m=l instabilities (ISaha et al.ll2007l) . In principle the lat- 
ter physical processes are likely to occur homogeneously for any 
sample independently of its environmental properties. Thus, the 
intrinsic asymmetry distribution found in our sample of isolated 
galaxies is likely due to a combination of these processes. Owing 
to the lack of spatial resolution in our data, at this moment we 



cannot distinguish their respective importance. High-resolution 
observations of isolated galaxies are a good probe to shed light 
into the origin of these more subtle asymmetries, though they are 
likely widespread in all kinds of environments. 

We have started a follow-up study of the origin of H I asym- 
metries using Very Large Array (VLA), expanded VLA and 
Giant Meter Radio Telescope (GMRT) aperture synthesis H I ob- 
servations of a subsample of ~ 20 isolated galaxies, which will 
be presented in a forthcoming paper. We selected galaxies cov- 
ering the wide range of asymmetries found in our sample. One 
of the isolated galaxies presenting an asymmetric profile in this 
subsample is CIG 96 (NGC 864), whose H I synthesis imagin g 
from the VLA has been studied in detail in lEspadaet al.l (l2005l) . 
The asymmetry in the H I profile is associated with a strong kine- 
matical perturbation in the gaseous disk of the galaxy, where on 
one side the decay of the rotation curve is faster than Keplerian. 
Although a companion is detected, no tidal tail is found, and it 
is probably not massive enough to have caused this perturbation. 
Probably we are witnessing the recent merger of a small gaseous 
companion. 



7. Summary and conclusions 

We used H I global velocity profiles for a large sample of isolated 
galaxies to /) quantify the rate and amplitude of H I asymmetries 
in isolated spiral galaxies, where environmental processes such 
as tidal interactions and ram pressure are minimized, //) study the 
role of the environment on the H I asymmetries, and ///) study 
possible correlations between HI lopsidedness and the prop- 
erties of the stellar component, including their morphological 
types, signatures of optical perturbation, bar and spiral strengths, 
as well as optical and FIR luminosities. 

To quantify the HI asymmetry, we calculated a flux ratio 
asymmetry parameter {Afiux ratio)- We restricted our study to a 
sample oi N = 166 galaxies (the HI refined subsample) for 
which we minimized undesired artificially induced lopsidedness 
by avoiding large uncertainties owing to the rms of the profile, 
determination of the mean velocity, and pointing off'sets. 

We found that a half-Gaussian curve properly fits the 
A flux ratio distribution of this refined sample, with a cr = 0.13. 
We suggest that if we deconvolve other sources of errors such as 
baseline fitting and random pointing off'sets, then the underlying 
cr is reduced to cr ~ 0.1 1. We confirm that by using this sample 
we eff'ectively minimize nurture eff'ects, because there is a lack 
of correlation between H I asymmetries and isolation parameters 
such as tidal force (one-on-one interactions) and number density. 

We compared the distribution of HI asymmetries of pre- 
viously studied field galaxies with that of our isolated galaxy 
sample. A half-Gaussian fit does not successfully reproduce in 
general the asymmetry distribution of field samples. Indeed, the 
intrinsic cr is larger in field samples than in isolated galaxies. 
This is likely a result of the lack of an isolation criterion in the 
selection of field galaxies, which are likely contaminated by in- 
teracting objects. This suggests that environmental mechanisms 
(producing short-lived eff'ects ~ 1 Gyr) are fundamental mech- 
anisms to produce H I asymmetries in galaxies, and are indeed 
responsible for the ~ 10-20% diff'erence in the HI asymmetry 
rate we see in field galaxies with respect to isolated galaxies. 
The asymmetry distribution of galaxies in denser environments 
is likely even wider and more skewed. 

Within the isolated galaxy sample, we did not find any strong 
correlation between the HI asymmetry and internal properties 
such as the morphological type, optical and FIR luminosities or 
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signature of interaction. A signature of perturbed optical emis- 
sion is not a necessary condition for the H I profile to be asym- 
metric, and vice versa. We found a trend for larger H I asymme- 
tries to be located in more FIR luminous galaxies that are likely 
interacting objects. We also found evidence that galaxies with 
higher spiral arm strength have lower HI asymmetries. 

The here presented HI refined subsample can be used to 
study the origin of intrinsic H I asymmetries in isolated galax- 
ies, and it is also a baseline for samples of galaxies in denser 
environments with H I data that are properly evaluated for in- 
strumental eff'ects. This can help to shed light into the relative 
importance of diff'erent environmental and internally generated 
processes in shaping the H I disks of galaxies. 
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Table 1. Asymmetry quantification for the H I sample 



CIG 


Visual classification 


^flux ratio 


AA(rm5) 


I^{mean vel.) 


^^-^flux ratio 


2 


1 


1.009 


0.065 


0.029 


0.071 


4 


1 


1.079 


0.026 


0.015 


0.030 


8 





1.080 


0.028 


0.018 


0.034 


9 





1.108 


0.027 


0.016 


0.032 


11 


2 


1.456 


0.058 


0.037 


0.069 



Note. 1) CIG number, 2) visual classification (0 = symmetric, 1= slightly asymmetric, 2 = asymmetric), 3) flux ratio asymmetry pa- 
rameter A flux ration 4) AA(rm^): uncertainty owing to the rms of the HI profile, 5) AA(mean vel.): Afiuxratio uncertainty owing to the 
determination of the mean velocity, and 6) AAfiux ration the final derived uncertainty of Afiuxratio^ including the effect of pointing off- 
sets. The full list is available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5), via 
http : //cdsweb . u-strasbg . fr/cgi-bin/qcat? J/A+A/vvv/ppp or from http : //amiga . iaa . es. 



Table 2. Comparison between the visual classification and A 



flux ratio 



Visual classification N Mean Median cr 

Symmetric 141 1.08 1.07 O.OT" 

Slightly asymmetric 126 1.13 1.13 0.09 

Asymmetric 45 1.37 1.32 0.26 



Table 3. Comparison of the half-Gaussian cr and H I asymmetry rate between samples of isolated/field galaxies 



Sample 


N 


cr 


Afluxratio> 1-26 


H I refined sub sample 


166 


0.13 


9% 


Haynes et al. (1998) 
Haynes et al. (1998) no CIGs 
Matthews et al. (1998) 
Bournaud et al. (2005) 


104 
80 
30 
76 


0.13 
0.13 

0.23 


9% 
10% 

17% 

22% 


All, no CIGs 


186 


0.17 


16% 
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Fig. 1. Basic properties (normalized distributions) of the H I sample (N = 31 2 galaxies, blue dashed line), the H I refined subsample 
(N = 166 galaxies, blue solid line), as well as the optically complete sample (IVerdes-Montenegro et al.ll2005l red dash-dotted line): 
a) velocity (V[ km s"^]), b) morphology (T(RC3), as in the RC3 catalog), c) optical luminosity (log(L5[Lo])) and d) FIR luminosity 
(log(L™[Lo])). 



GIG 266. A = 1.05±0.04 




6300 6400 6500 6600 6700 6800 6900 7000 
Velocity{km s^^) 



Fig. 2. Example of a symmetric H I profile: CIG 266, Afiux ratio = 1 .05 ± 0.05. The points where the horizontal (black) line intersects 
the profile correspond to the low (V/) and high (Vh) velocity ends at a 20% level with respect to the peak. The derived mean velocity 
(Vm) at a 20% level is plotted as a (red) point. 
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GIG 421. A = 1.15±0.02 




6300 6400 6500 6600 6700 6800 6900 7000 
Velocity{km s^^) 



Fig. 3. Example of a slightly asymmetric H I profile: CIG 421, Afiux ratio = 1.15 ± 0.03. See description in the caption of EigureO 



GIG 361. A = 1.51 ±0.03 




7000 7200 

Velocity{kva s^^) 



Fig. 4. Example of a strongly asymmetric HI profile: CIG 361, A fiux ratio = 1-51 ± 0.03. See description in the caption of Eigure|2l 
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Fig. 5. Upper panel) The Afiux ratio distribution (solid line histogram) of the HI sample (N = 312) and its best half-Gaussian fit 
(dashed line). Lower panel) The residual of the half-Gaussian fit to the observed Afiux ratio distribution. 
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Fig. 6. Upper right) Uncertainty distribution owing to the rms of the HI profiles, AA(rms). Upper left) Uncertainty distribution 
produced by errors in the mean velocity, AA(mean vel). Bottom) Uncertainty distribution of A fiux ratio for the HI sample (N = 312) 
combining the eff'ect of AA(rm^), AA(mean vel) and the small contribution of AA(pointmg offset) (See Sect. 13. 2. it . The H I refined 
subsample (AAfiux ratio < 0.05, indicated as a dotted line in the lower panel) is shown in the plots as blue filled histograms. Best 
Gaussian fits are presented for all distributions as dashed lines: for the HI sample: a) yu=0.02, (T=0.02, b) ju=0.04, (T=0.02, and c) 
yu=0.04, o-=0.03; and for the H I refined subsample: a)ju=0m2, cr=0.007, /?)yu=0.029, o-=0.012, and c)ju=0.033, o-=0.012. 
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Fig. 7. Comparison between the visual classification and the Afiux ratio parameter for the N = 312 galaxies in the H I sample (sym- 
metric: white histogram, slightly asymmetric: blue or dark gray histogram, and strongly asymmetric: green or light gray histogram). 
Mean and standard deviations (cr) for each distribution are also shown by vertical and horizontal lines, respectively. 
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Fig. 8. A flux ratio distribution for the diff'erent cuts in AAfiux ratio from 0.03 to 0.1 1 in bins of 0.02. Note that AAfiux ratio < 0-05 (N 
= 166, purple filled histogram) corresponds to the HI refined sample (see Sect. 13.41) , and the HI sample (N = 312) is the blue solid 
line histogram. 
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Fig. 9. Upper panel) Afiux ratio distribution of the HI refined subsample (AAfiux ratio < 0.05) (N = 166, blue filled histogram), in 
comparison with that of the HI sample (N = 312, solid line histogram) (see Sect. 13.4b . A half-Gaussian fit (red dashed line) to the 
HI refined subsample is presented. The half-Gaussian curve is characterized by a standard deviation cr=0.l3, lower panel) The 
residual of the half-Gaussian fit to the observed Apux ratio distribution for the H I refined subsample. 
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Fig. 10. a) Distribution of inclinations for the N = 3 12 galaxies in the H I sample (solid line histogram) and the H I refined subsample 
(filled histogram), b) Afiux ratio versus inclination, from / = 10 to 90° in 10° bins for the HI refined subsample. Red points and their 
error-bars indicate the mean (connected by red solid line) and standard deviation, and green points the median and the median 
absolute deviation. 
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1.2 1.3 

^flux ratio 



Fig. 11. Isolation parameters vs Afiuxratio- The isolation parameters (IVerley et al.ll2007bb are the local number surface density 
parameter tjk to the K-th neighbor, where K = 5 (upper panel), and the tidal s trength parameter Q (bottom panel), which only takes 
into account similarly size neighbors (factor 4 in size, as defined in iKarachentseva. 1973.) . 
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Fig. 12. Compa r ison o f the normalized Afiuxratio distribution between our HI refined sample (black solid histogram) and 1 ) 
l2005l) (re d dotted-dash e d hist ogram), 2) a combined sample including HI data in iBournaud et alJ (l2005l) . 



and iHaynes et alJ (1 19981) excluding CIG galaxies (green dotted line). Solid curves are the half-Gaussian 



Bournaud et al. 



Matthews et al 



curves fitted to each distribution. See Table[3]for a comparison of the cr's of each half-Gaussian curve. 



D. Espada et al.: The AMIGA sample of isolated galaxies. VIIL t 



17 



^ 1.0 
g 0.8 
^0.6 

■I 0-4 
^ 0.2 



on 





1 1 1 1 1 1 1 

a) _ 


- 


Matthews et a/. (1998) _ 

1 —This work 


'-^. _ _ 





1.0 
0.8 
0.6 
0.4 
0.2 



,0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1. 



O.Q 



L 


1 1 1 1 1 1 1 

b) _ 

-1 

■■ 1 Bournaud et al.{2005) _ 


*~| ~ -1 



1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.^ 



;5 




- Haynes et aL( 1998) 
{no CIGs) 



1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1. 

^flux ratio 



1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.^ 

^flux ratio 



Fig. 13. Comparison of th e cumulative Afiux ratio distribution between our samp le (black solid line) a nd other samples (red dashed 
lines): g) .Matthews etalJ (1998 , A^ = 30), /?) BoumaudetaL ( 2005, A^ = 76), c) lHavnes et ^1(119981 A^ = 104)) (A^ = 106), and d) 
iHaynes et alJ (Il998l) excluding CIG galaxies {N = 80). See Table[3]for an asymmetry rate comparison at a Apux ratio = 1-26 level. 




T{RC3) 



Fig. 14. A flux ratio and T(RC3) (ISulentic et alJ 120061) for the HI refined subsample. Red points and their error-bars indicate the 
mean (connected by red solid line) and standard deviation, and green points the median and the median absolute deviation, for each 
morphological type from T(RC3) = -5 to 10 (E to Im). 
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Fig. 15. Left) A flux ratio versus log(Lfi[Lo]). Symbols are as in Figure [141 Right) Cumulative probability distribution of A fiux ratio 
for log(Lfi[L0]) < 10 (solid line) and log(Lfi[Lo]) > 10 (dashed red line), using the HI refined subsample. 
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Fig. 16. A flux r atio ve rsus spiral strength (Qs) for 40 CIG galaxies overlapping between the HI sample and the GIG galaxies in 
iDurbala et alJ (l2009l) . The fit to the data points (slope and intercept are -0.34 and 1.16, Pearson's correlation coefficient p = -0.45) 
is shown as a (red) solid line. The six outliers have been ignored in this fit. 
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Fig. 17. Left) A flux ratio versus log(L/7//?[Lo]). Symbols are as in Figure [T4l Right) Afiux ratio cumulative probability distribution for 
three Lfir bins: 9 < log(L™[L0]) < 9.5 (solid line), 9.5 < log(L™ [L©]) < 10.0 (blue dashed line) and 10.0 < log(LFiR[LQ]) < 11.0 
(dash-dotted line). The (green) dotted line represents the Afiux ratio cumulative distribution for those galaxies with an upper limit. 



